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THE  UNDERPRESSURE  FIELD  FROM  EXPLOSIONS  IN  WATER 
AS  MODIFIED  BY  CAVITATION 


1.  INTRODUCTION 

In  connection  with  a  current  study  of  the  effects  of  underwater  explosion 
pressure  waves  on  marine  organisms,  it  was  desired  to  determine  the  effect  of  bulk 
cavitation  on  the  pressure  field.  In  particular,  since  the  negative-going, 
surface-reflected  shock  wave  has  been  suggested  as  an  important  element  in  damage 
production  to  some  classes  of  fish,  the  prediction  of  these  negative  pressures  was 
undertaken.  Existing  cavitation  theory  has  been  used  to  map  the  cavitated  region 
for  several  geometries  of  interest.  Negative  pressure  amplitudes  beneath  the 
cavitated  region  have  been  calculated  using  two  different  approachs  in  addition  to 
that  suggested  by  Arons,  et  al,  (Ref  1);  the  results  of  the  three  approaches 
are  compared  with  some  existing  data. 

2.  CAVITATION  THEORY 

If  cavitation  is  not  considered,  the  amplitude  of  the  surface  reflected 
shock  wave  is  usually  calculated  by  assuming  a  negative-going  image  source  of 
strength  equal  to  the  actual  source.  Since  the  occurence  of  cavitation  markedly 
reduces  the  strength  of  the  reflected  wave,  non-cavitation  predictions  are  unrealistic. 

The  history  of  cavitation  theory  shows  no  lack  of  effort,  especially  since 
World  Nar  II.  References  1  through  8  are  only  a  few  of  the  papers  existing  which 
cover  aspects  of  cavitation  theories.  In  actuality,  the  groundwork  for  all  these 
theories  was  laid  by  Kennard  in  two  papers  in  19^3  While  there  have  been 

refinements  and  additions  to  this  work  (especially  as  concerns  the  dynamics  of  the 
cavitated  region  after  initial  formation,  which  is  of  little  interest  for  the 
application  of  this  report),  the  basic  assumptions  bearing  on  the  formation  of  the 
cavitated  regi  m  have  remained  unchanged  (as  well  as  untested  by  experiment). 

The  development  presented  below  follows  that  of  Arons,  et  al  This  t 'eatment  was 
chosen  for  its  conceptual  and  computational  simplicity. 

Reference  1  makes  three  basic  assumptions: 

(a)  "Cavitation  begins  when  the  pressure  in  the  water  drops  to  a  fixed  value, 

Pc,  the  cavitation  pressure.  We  will  take  Pc  =  0  psi." 

(b)  "The  pressure  at  a  point  on  the  front,  of  the  reflected  wave  can  be 
described  as  emanating  from  an  image  charge  above  the  surface,  having  a  variable 
apparent  weight,  Wi ,  so  adjusted  as  to  account  for  the  pressure  dimhnvtion  due  to 
cavitation. " 

( c )  "Wi  is  a  continuous  non-increasing  function  of  distance  along  a  ray  away 
from  the  image  charge,  being  equal  to  the  true  charge  weight,  W,  before  any 
cavitation  takes  place,  decreasing  monotonically  in  regions  of  cavitation,  and 
remaining  constant  in  regions  of  no  cavitation." 
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The  total  pressure  in  the  water,  before  cavitation  has  occurred,  is  given  by 
the  sum  of  the  direct,  surf ace-reflected ,  and  hydrostatic  pressures.  Assuming  for 
calculations!  convenience  that  the  shock  wave  is  an  instantaneous  rise  to  a  peak 
pressure,  followed  bv  an  exponential  decay  of  time  constant  0,  the  pressure  at  a 
point  in  the  water  at  the  time  of  arrival  of  the  surface  reflection  is: 

P  e"t76  +  P  +  P  ,  (1) 

o  r  z 

where  Pq  is  the  shock  wave  peak  pressure,  t  is  the  time  delay  bet/een  the  direct 

and  surface-reflected  arrivals,  P  is  the  surface-reflected  peak  pressure,  and  P 

1*  z 

is  the  hydrostatic  pressure.  Actually  the  relationship  for  pressure  as  a  function 

of  time  in  the  shock  wave  is  only  goo!  until  the  time  is  about  one  0;  at  later 

times  the  pressure  may  be  well  above  the  calculated  value.  For  calculating 

pressures  in  the  immediate  vicinity  of  the  cavitation  region  it  is  probably  quite 

adequate.  According  to  assumption  (a)  we  have: 

P  e't70  +  P  +  P  =  0  (2) 

o  r  z 


The  coordinate  system  to  be  used  in  the  following  calculation  is  shown  in  Figure  1. 


Assuming  that  the  shock  wave  peak  pressure  may  be  calculated  from  a 
similitude  equation  of  the  form: 

yl/3  1,13 

P  =  K  (JS-jp)  , 


(3) 


we  have 


w173 

P  =  K  (-V-) 

o  Ri 


1.13 


and 


(>0 


,Wil/3 
P  =  K  (-i - ) 

2 


1.13 


(5) 


where  K  is  an  empirical  constant  depending  on  the  explosive  composition,  R,  is  the 
slant  range  from  charge  to  gage,  and  Rg  is  the  slant  range  from  image  to  gage.  In 

addition,  the  hydrostatic  pressure,  P^,  may  be  assumed  given  by: 

P  =  P  +  H  (R  cos  a  -  d),  (6) 

Z  8,  JL 


where  P  is  the  atmospheric  pressure,  d  is  the  charge  depth,  H  is  c  constant,  and 
8. 

a  is  the  angle  between  the  line  from  image  t,o  gage  and  a  vertical  line  through 
image  and  charge. 
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FIG.  1  COORDINATE  SYSTEM  FOR  CAVITATION  CALCULATION 
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Combining  equations  2,  k,  5,  and  6,  we  have: 


1.13  -t/e 


-  K  (~— )  +  P  +  H  (R  cos  a  -  d)  =  0 

R2  a  2 


.1/3 


1.13 


Solving  for  the  variable  image  weight, 
1.13 


i  Rr 


W.  = 
i 


K 


T  j ^ 

P  +  H  (R  cos  a  -  d)  +  K 

Q,  C  A 


1.13  ~t/G- 


3.00 


1.13 


(7) 


(8) 


The  arrival  difference,  t,  is  given  by: 


t  = 


R1  ~  R2 
C 


(9) 


where  C  is  the  velocity  of  sound  in  water.  The  cavitation  region  may  now  be 
mapped  using  assumptions  (b)  and  (c).  At  a  fixed  a,  allowing  Rg  to  vary,  the 

upper  surface  of  the  cavitated  region  is  found  by  setting: 

W  =  W  (10) 

The  bottom  of  the  cavitated  region  is  found  when 

dW. 

5C=0  (11) 


Following  rays  from  the  image  charge  of  different  angles  a,  the  entire  cavitation 
zone  may  be  mapped  out.  To  calculate  the  negative  underpressure  at  points  below 
the  cavitated  region  using  the  approach  suggested  by  Arons,  et  al,  the  image  weight 
corresponding  to  equation  11  is  used,  and  spherical  spreading  is  utilized  to 
extrapolate  from  the  underpressure  at  the  bottom  of  the  cavitated  zone  to  the 
desired  point  in  the  water. 

One  modified  approach,  Gaspin-Price  plane  wave,  assumes  that  the  amplitude  of 
the  negative  wave  is  determined  at  the  bottom  of  the  cavitation  zone  in  the  manner 
outlined  above.  The  wave  then  continues  to  propagate  radially  from  the  image 
source  but  as  a  plane  wave  with  its  origin  the  bottom  surface  of  the  cavitated 
region.  It  is  obvious  that  this  is  not  completely  realistic,  but  it  seems  to  be 
closer  to  reality  than  assuming  a  point  source  with  spherical  spreading. 
Specifically,  it  should  be  better  at  positions  beneath  the  central  portion  of  the 
cavitated  region  (such  as  f  in  Figure  2(a)).  The  assumption  of  plane  wave 
propagation  probably  becomes  less  valid  when  the  radial  line  from  the  image  source 
passes  through  the  extremities  of  the  cavitation  region  (as  at  g  in  Figure  2(a)). 

The  duration  of  the  negative  phase  is  determined  by  the  time  that  the 
cavitation  remains  open,  i.e.,  until  cavitation  closure  occurs  along  the  ray 
connecting  the  image  to  the  gage. 
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There  are  undoubtedly  limits  to  the  regions  below,  and  just  out-from-under , 
the  cavitation  zone  in  which  the  modified  approach  may  be  used.  The  limitations 
have  to  be  set  by  appropriate  analytical  and  experimental  methods. 

Another  approach,  the  hydrostatic  origin  approximation,  assumes  that  the 
calculated  negative  pressure  at  the  bottom  cf  the  cavitation  region  directly  above 
the  point  of  interest  (a  gage  for  example)  is  propagated  downward  as  a  plane  wave. 
Essentially  the  pressure  due  to  water  and  atmosphere  above  the  cavitation  region 
is  removed  from  the  water  column  below  the  cavitation.  The  pressure  is  removed 
until  the  loading  pressures  and  cavitation  closure  restore  conditions  to  a  more 
nearly  normal  state.  The  long  duration  and  even  pressure  in  the  early  parts  of 
the  negative  phase  are  thus  explained. 

This  approximation  must  break  down  at  great  depths  and  at  positions 
horizontally  beyond  the  cavitated  region. 

3.  MODIFICATION  OF  THE  SURFACE  REFLECTED  WAVE  FORM  RY  CAVITATION 

Figure  2  illustrates  the  manner  in  which  the  wave  form  of  the  surface 
reflected  shock  wave  is  modified  by  cavitation.  The  five  pressure-time  curves  are 
lettered  to  correspond  to  the  five  successively  deeper  points  indicated  in  the 
sketch.  At  point  "a",  the  surface  reflection  does  not  lower  the  pressure  as  far 
as  absolute  zero,  and  no  cavitation  occurs.  At  this  point,  the  surface  reflected 
pressures  may  be  reasonably  calculated  assuming  an  image  source  of  equal  strength 
to  the  actual  source,  and  a  negative  exponential  wave  form.  At  point  "b",  on  the 
upper  edge  of  the  cavitated  zone,  the  reflected  pressure  is  just  sufficient  to 
lower  the  total  pressure  to  absolute  zero  and  cavitation  occurs.  After  the 
pressure  reaches  zero,  it  stays  relatively  stable  at  that  level  and  eventually 
merges  into  succeeding  events  (closure  ).  At  points  "c"  and  "d",  within  and 
on  the  lover  boundary  of  the  region,  the  reflected  wave  form  shows  a  sharp 
negative  spike  below  absolute  zero.  This  net  tension  results  from  the  finite  time 
necessary  for  the  water  to  cavitate.  On  many  recordings,  this  spike  is  not 
apparent  due  to  inadequate  frequency  response.  Following  this  spike,  the  pressure 
returns  to  a  more  or  less  stable  "plateau"  pressure  of  absolute  zero.  As  the 
wave  propagates  to  point  "e".  the  wave  form  retains  the  general  features  which 
characterized  "c"  and  "d",  but  the  reflected  pressure  is  no  longer  adequate  to 
produce  cavitation. 

The  underpressure  at  points  b,  c,  and  d  may  easily  be  calculated.  The  total 
pressure  at  these  points  reaches  absolute  zero,  so  the  underpressure,  is  just 
equal  to  the  hydrostatic  pressure  at  these  depths.  At  point  "e",  the  pressure 
may  be  calculated  either  using  the  image  weight  obtaining  at  point,  "d",  and 
allowing  for  diminution  of  the  underpressure  due  to  spherical  spreading  between 
points  "d"  and  "e",  or  by  calculating  the  underpressure  at  "d"  and  assuming 
that  this  underpressure  propagates  as  in  plane  wave  to  "e."  We  have  not  found 
data  to  check  conditions  at  a  point  such  as  "e"  which  is  not  under  the  cavitated 
region. 

U .  RESULTS 

Computer  programs  have  been  written  to  calculate  the  negative  pressure  below 
the  cavitated  region  and  the  boundary  of  the  region;  they  are  listed  in  the 
appendices.  They  are  written  in  BASIC  and  have  been  run  on  the  CDC  twOO  computer. 
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The  cavitated  region  has  been  mapped  for  several  conditions  of  interest,  as 
sacra  in  Figures  3  through  6.  As  mentioned  above,  this  calculation  is  essentially 
unchecked  by  experiment-  The  experiment  described  in  reference  6  roughly 
confirmed  the  calculation  of  the  upper  boundary  of  the  cavitated  region,  but  to 
the  authors  knowledge,  the  calculation  of  the  lower  boundary  has  never  been 
experimentally  verified.  Ho  data  is  presently  available  to  make  a  direct 
verification. 


Data  is  available,  however,  to  compare  with  calculations  of  the  magnitude  of 
the  underpressure  below  the  cavitated  region.  Computer  programs  were  written  to 
make  this  calculation,  and  several  samples  of  available  data  were  analyzed  to 
determine  the  negative  "plateau"  underpressure.  'These  measurements  are  of  a 
rather  crude  nature,  since  the  instrumentation  was  set  up  to  record  the  much 
higher  amplitude  of  the  direct  shock  wave,  and  the  negative  pressures  are  a  small 
fraction  of  the  full  scale  reading.  They  do,  however,  give  some  indication  that 
the  calculations  of  underpressures  are  valid,  particularly  if  plane  wave 
propagation  is  assumed,  and  thus  provide  an  indirect  verification  of  the 
calculation  of  the  lower  cavitation  boundary.  Table  I  and  Figures  7  through  10 
present  this  comparison.  Also  indicated  in  Table  I  are  the  corresponding  under¬ 
pressure  as  calculated  by  assuming  an  image  source  of  strength  equal  to  the 
actual  charge  (no  provision  being  made  for  cavitation).  It  is  clear  that  the 
cavitation  calculations  give  a  much  more  accurate  representation  of  the  measured 
pressures.  Sir-e  data  from  charges  of  0.U9  and  It  x  10^  pounds  are  reasonably  well 
predicted,  ti  present  calculation  is  considered  adequate  for  purposes  of 
experimental  design.  A  more  extensive  comparison  would  be  desirable,  however. 

The  fact  that  these  pressures  are  reasonably  well  modeled,  and  that  the  position 
of  the  lower  boundary'  influences  these  pressures,  indicates,  albeit  inuirectly, 
that  trvre  probably  is  some  validity  to  the  lower  boundary  calculation. 

5.  CONCLUSION 

Existing  bulk  cavitation  ^heory  may  be  used  to  predict  the  extent  of  the 
cavitated  region  caused  by  underwater  explosions  and  also  the  negative  pressure 
amplitudes  below  tne  cavitated  region.  While  the  calculation  of  the  extent  of  the 
region  is  largely  untested  by  experiment,  preliminary  comparisons  with  data 
indicate  that  these  calculations  of  the  underpressure  amplitudes  are  a  considerable 
improvement  over  calculations  made  without  considering  cavitation. 
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FIG.  5  BOUNDARY  OF  CAVITATED  REGION 
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FIG.  6  BOUNDARY  OF  CAVITATED  REGION 
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TABLE  I.  UNDERPRESSURE  COMPARISONS  (continued) 
Pressures  in  psi 


O  CO  1A  H  KO 

iA-3-_=tiA.=r.j-oo*  *  * 

i  t  i  t  i  i  i 


0\  H  t- 
H  CM  H 

I  I  I 


\OCO  Ol  t“4-  W  J-  iAt-0\ 

\o  m  hcov)  j  w  o  o\co 

-3-  4-.3-OOOOOOOOOOCMCM 

I  I  I  I  I  I  I  I  I  I 


-=r  oo  oo 

-3-  oo  co 

CM  H  I 

I  I 


i  r-J  CM 

O  IT 

i  --H  i 


CM  0\  t-MO  tA  00  OO  CM  CM  H 


I  I  I  I 


cr\  CO 

rH  rH 
I  I 


C0t-00O\C^lA^J-CMCMH  tA  00 
00  OJOJHHHHHHH  H  i— I 

i  I  I  I  I  I  I  I  1  I  I  I 


cooooocooocooocooooo 

i  i  i  i  i  i  i  i  i  i 


o\  oo 

rH  rH 
I  I 


O  H  CM  00.3  IA  IA  C-COOO 
CMC0_3OV0CM00-3l'-O 
CM  CM  00-3-3  IA  IAVO  VO  C~- 


a 

bO 

W 

o 

C 

-3 

CO 

-3- 

o 

•H 

•rH 

$3 

oo 

-=f 

CO 

o 

-p 

O 

ir\ 

1 — 1 

ctf 

a) 

•rl 

6 

<l> 

•p 

•H 

Sh 

cd 

x 

ft 

•P 

0 

CO 

•rl 

P 

0) 

> 

ft 

rH 

> 

a) 

ft 

a) 

a 

O 

< 

U 

> 

•H 

bO 

a 

u 

0) 

a 

o 

o 

O 

o 

•rl 

<v 

a 

♦rl 

o 

rH 

rH 

rH 

bO 

33 

<3 

u 

CM 

•rH 

ft 

rH 

0) 

!U 

CO 

ft 

Tj 

O 

•rH 

C 

a) 

(/) 

O 

•H 

V 

S3 

•rH 

ft 

•rH 

0 

P 

w 

u 

a} 

0} 

ft 

■P 

O 

1 

•p 

VO 

LA 

o\ 

w 

1 

c 

3 

o 

O 

-3 

0 

w 

•rH 

o 

• 

• 

• 

c 

ft 

33 

O 

iH 

o 

o 

w 

■P 

O 

a) 

•h 

O 

33 

<3 

O 

o 

H 

«■—- 

/ — ' 

«- — > 

ctf 

& 

O 

' — ' 

' - ' 

13 


Distorted  by  the  Bottom  Reflection 
No  Data  or  Not  Calculated 


NEGATIVE  PRESSURE,  PSI,  RELATIVE  TO  AMBIENT 
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FIG,  7  MEASURED  VS  CALCULATED  UNDERPRESSURES 
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00100 
00110 
00120 
00130 
001  40 
00150 
00160 
00170 
00180 
00190 
00200 
00210 
00220 
00230 
00240 
00250 
00260 
00270 
00230 
00290 
00300 
00310 
00320 
00330 
00340 
00350 
00360 
00370 
00380 
00390 
00400 
00410 
004?0 
00430 
00440 
00450 
00460 
00470 
00480 
00490 
005  00 
005  10 
00520 
005  30 
005  40 
00550 
00560 


APPENDIX  A 

COMPUTER  PROGRAMS  FOR  CALCULATING.  BOUNDARY  OF  CAVITATED  REGION 


REM  READ  K,  ALPHA  F3  R  PRESSURE 
READ  K*  A2 

REM  READ  K,  ALPHA  FOR  THETA 
RFAD  A3,A4 

REM  ARONS  METHOD  FOR  CAVITATION 
REM  W=  CH G  WT,  LBS. 

REM  D=  CHG  DEPTH,  FT. 

RFM  R8  =  H3  RI  Z  3N  TAL  RANGE  INCREMENT 

REM  R9=MAX  H3KIZ3NTAL  RANGE 

READ  W,  D,  R8,  R9 

LET  Pl=  14.7 

LET  H=.  45 

LET  C=  5000 

PRl  N  T"W=  "W,  "DE  PTH=  "D 

PRINT 

REM  PRINTOUT  IS  HO R.  RANGE  AND  DEPTH  OF  CAV  START  AND  STOP 
REM  ALONG  A  RAY  FROM  IMAGE,  AND  MIN.  OF  IMAGE  WEIGHT 
PRINTER  START",  "Z  START",  "R  STOP",  "Z  STOP",  "IMAGE  WT" 

PRINT 

W2=Wt  (A2/  3) 


H  1  =  -  R8 
1  =  0 

H  1  =  H  1  +  R8 
A=pTNCH1/D) 

IF  H  1  >  R9  THEN  01070 
LET  R1=D/C0S(A) 

REM  FIND  W<  IMAGE)  =K 
L  FT  1  =  1+1 
GO  SUB  0063  0 


Reproduced 

best  available_coHiL_ 


RFM  '  X  =  WC  IMAGE) 


LFT  F=X-W2 

IF  AB  S(  F)  <  1  F-  6  THEN  00480 
IF  I  >  5  0  THFN  00460 
GO  SI 'R  00730 
LFT  R 1 =R1 - (F/X l ) 

GO  TO  00370 

PRI N T"NO  CONVERGENCE  ON  CAV  START" 

GO  TO  01070 

LFT  0 1  =  R 1  *  CO  SC  A  )  -  D 

LET  G2=  Rl  *  SI  N(  A) 

RFM  FIND  MIN  OF  WCIMAGE) 

LET  1  =  0 
LET  1=1+1 
GO  SUB  0078  0 

IF  ABSCX1XIF-6  THEN  00610 
IF  I >5  0  THEN  00590 
GO  Si  'P  008  60 


A-l 
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APPENDIX  A  (continued) 


00570 
00580 
00590 
006  00 
006  10 
00620 
006  30 

006  40 

006  5  0 
00660 
00670 
00680 
00690 
00700 
007  10 
00720 
007  30 
007  40 
007  50 
00760 
00770 
00780 
00790 
008  00 
008  10 
008  20 
008  30 
008  40 
00850 
00860 
008  7  0 
00880 
00890 
009  00 
009  10 
009  20 
009  30 
009  40 
009  50 
00960 
00970 
0098  0 
00990 
01000 
01010 
01020 
01030 
01040 
01C50 
01060 
01070 


LET  H 1  =  R 1  -  C  X  1/X2) 

G3  T0  005  20 

PRIN  T"N0  C3NV  3N  CAV  ST3  P" 

G3  T0  01070 

LFT  03=  R1  *C3  SC  A)  -  D 

LET  G4=R1*SINCA) 

Q5=Xt  C  3/  A 2) 

PRINT  Q2»  Ql>  Q4>  G3»  Q5 
PRINT 

G3  T3  00310 

REN  CALCULATE  WC  IMAGE) 

LFT  R2=  FN  AC  R 1 ) 

LFT  B=FN PC  R2) 

LFT  T=  C  R 1  -  R2) / C 
LFT  G  1  =  -T/R 
LFT  G2=FXPCGl) 

G3=G2*X*W2/<R2»  A?) 

LET  G4=H*C  R1*C3  SC  A)-D) 

X=CRlt  A2)*C  Pl+G4+G3)/X 
RETURN 

REM  CALCULATE  DERI  V  0F  WC IM AGE) . . • X  1 
G0  SUB  0068  0 
LFT  Dl  =  FN CC  R 1 ) 

LFT  D?=  FN  DC  R 1 ) 

L  1=  C  C  R2/B)*D1)+A2*D2 
LFT  L2=L1*G3/R2 

X  1  =  C  C  R 1  t  A2>/K)+C  m<C3SC  A) -L 2)  +  C  A2*X/ R 1  > 

rfturn 

RF.M  CALC 'LATE  2ND*  DERIV  OF  WC  l M  AGE) ...  X 2 
D4=FNFC  Rl) 

LFT  D 3=  -  D 4/  C 

Ml  =  -C  R  It  A2)  *W2*G?/C  R2»  C  1+A2>) 

*2=  C  R2*D3/B)  +  C  C  D 1  *  D2)  /  B)  +  A2*  D4 

M3  =  C  F2*  D 1/  B)  A2*  D2 

M  4=  C  C  R2*D1/P)  +  0C  A2+  1)  *D2)  )/R2 

M  5=  A2*  R 1 1  C  A 2-  1 )  *C  M*C3  SC  A)  -  C  G3/F2)  3) /* 

LET  M6=  C R1*X 1-X)/ C Rlt 2) 

M6  =8  6*  A? 

LFT  X2=M  1*C>M2-CM3*M4)  )+M5+M6 
Rl TURN 


RFM 

VN  A=  R2 

DFF 

FNACR1) 

=  SQRC  R 

1  r  2+  4*  D*  2-  4*  D* 

R 1  *  03  SC  A) 

REM 

FN  P=  THE 

TA 

DFF 

FN  PC  R)  = 

A3* 1E- 

3*W 

2*C  W2/R)t  A 4 

i 

REM 

FN  S  CL  D 

,  E  ARF 

IN 

TFRMFDI  ATF 

CALC 

S 

DFF 

FNCC  RI) 

=  c  t-c c 

Rl- 

2*D*  DISC  A)  ) 

/  R2) 

)/C 

OFF 

FN  DC  R  1 ) 

=  C  Pl-2 

*  D* 

C3  SC  A)  )  /  R2 

DFF 

FNFCR1) 

=  C  R2-C 

Rl- 

2*  D*  C0  SC  A) ) 

*  D2) 

/R2t  2 

DATA 

.  20800, 

i  •  1  3*  • 

06, 

18 

DATA 

4C0C0, 

POO,  50 

,  1  1 

OC 

END 

j2  0 

A-2 


HOLTR  72-103 


APPENDIX  B 

COMP OTffl  PROGRAM  FOB  CALCULATING  UNDER  PRESSURE  BELOW  CAVITATED  REGION 

C010C'  FFM  GASPIN- PRICE  PLANF  WAVE  PROPAGATION 
C01IC  PEH  ALONG  RAYS  FROM  IMAPF  SIPRCF. 

00120  REM  KFAD  ALPHA  FT P  PRFSSURF 

00130  FFH  MSF  CONSTANTS  TO  ORIGIN  PRESSURE  IN'  PSI  • 

001 AO  RFAD  K*A2 

C01S0  PFH  RFA1?  X#  ALPHA  FOR  THF.TA 

CGI  AO  REM  "SF  CONSTANTS  TO  1PTAIN  TIME.  CONSTANT  IN  MILLI  SECONDS* 
0C170  RFAD  A3*AA 

00130  REM  PROGRAM  CALCILATFS  UNDERPRESSURE  F0R  GAGE  L  T  CA  TI  ON 
C0190  PFM  HELOP  CAVITATED  ZONE 

00200  RF.M  AFDNS  METHOD  FTP  CAVITATION 
00210  RFM  V:=CHG  WT>  LGS. 

C0220  PF-*  D=CHO  DEPTH#  FT*  • 

00230  REN  H2=PFS1RF0  GAG17  HORIZ  RANGF 
002*0  REM  09=  DF SI  RED  GAGE  D2PTH*ET. 

00250  RFAD  V’»  D 
C02A0  Pnl N T"W=  "K#  "DFPTH=  "D 
00270  RFAD  H2,  D9 
00230  LFT  Pl=  1-6.  7 
00290  LET  H=.  -65 
0030C  LFT  C=5000 
00310  PRINT 
CG320  '/  p=  Vt  (  A 2/  3) 

00330  A=ATNCH2/<D+D9)) 

003*0  LEI  R1=D/C*S<A> 

00350  1=0 
00360  LFI  1=1+1 
00370  GOSMR  0035  0 
0033 C  LFT  F=X-W2 
00  390  IF  AHSCFX  1F-6  THFN  00A30 
00-600  IF  I>50  THEN  00  A  AO 
OOAIC  FT  SOP  009  AO 
00A20  LET  R 1=  P.l-  C  F/X  1 ) 

00A30  GO  T-*  00360 

OOA^O  PKINT"NO  CONVEHC-FNCE  ON  CAV  START" 

00*50  GO  TO  09  000 

O0A60  PRINT  "GAGE  APOVF  CAVITATFD  RFGI0N"JG1 
00A70  GO  TO  005  00 
00*30  LFT  0 1  =  R1  + CT  S(  A) -  D 
C0A9  0  IF  D9<  6 1  THEN  00A60 
00500  LFT  02=  PI*  SI NC  A) 

005  10  LFT  1=0 
00520  LFT  1=1+1 
005  30  GTS'lR  009  AO 
005  AO  IF  APSCX1X1F-6  THEN  00630 
00550  IF  I > 5 0  THFN  00590 
00560  GO  SUP  01010 
00570  LFT  P 1  =  K 1  -  C  X  1  / X 2) 

005 HO  GO  TO  005  20  9  t 


B-l 
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APPENDIX  B  (continued) 


0059  0 
00600 
00610 
C0620 
00630 
006  40 
00650 
00660 
00670 
006  8  0 
0069  0 
007  00 
00710 
00720 
007  30 
007  40 
0075  0 
0076  0 
00770 
0073  0 
0079  0 
008  00 
008  1  0 
003  20 
003  30 
008  40 
00350 
00360 
0037  0 
0038  0 
0089  0 
009  00 
009  10 
0C9  20 
009  30 
009  40 
009  60 
009  60 
009  7  0 
00930 
0099  0 
01000 
01010 
01020 
0103C 
0  1 0  40 
01050 
0 1 06  0 
01070 


PRIMT"\n  C3NV  3M  CAV  STOP" 

G3  TO  09000 

PRINT  "GAGE  IN  OR  ARTVF  OAVITATFP  REGION"? 03 
GO  T3  00650 
LET  03*  R 1  *  CM  S<  A)  -  R 
IF  D9<  03  THEN  00610 
LFT  04=  Rl*  SIN  (A) 

05=  X  t  (  3/ A  2) 

PRINT  "IMAGE  CHARGE  V*TGHT"05 
REM  CALCULATE  » 'N  EPRE  S  SI  IRE 
R  3=  h  1 

Y  1  =  K*X/ <  R 1 1 A2) 

M7=D-03 
R4=C 
T=CR 
M8=  - 
M9=  FX  PC  M3 ) 

R=,vi9  *!<*'.•  2/C  R4t  A2) 

P9=P-Y  1 

I  r*  CAlCI'LATF  UNDrF  Pt.FSSURF  HT'HIUT  CAV/I  T&TI  IN 

Y  2=8 2/  C  1-3*  A?) 

P8=  P-  Y2 

PRINT’T  GAGE="U9j  "H  OAOF="HP>"P  N  FG=  "P9,  "PN  EG  NMCAV=',P8 

Peel  NT"  P<  T)  =  "RJ  "T/THFTA="M8 

PRINT 

CO  T")  00270 
LET  R2=  FN  AC  r  1 ) 

L FT  -=FNPCF'?) 

LFT  T=<  rtl-h«)/C 
L  FT  0  1  =  -  TV  R 
LET  G?=  rA  P'  G  1 ) 

G3=  02  t=V'2/  (  RP»  A2) 

L  ~T  G 4=  H*  (  K  1  •* CO  S(  A  >  -  0) 

X=C  Rlr  A?)  K(  Fl  +  G4+G3)/-<  • 

RFTMRN 
G 1  SUE  08350 
LET  niaFNC(hl) 

LET  '  D2=  FN  nC  »<!) 

L  1  =  CC  R2/R)  M)1>  +  A2*R? 

LFT  L  2=L  1  *03/ H2 

X  1  =  (  (  R 1 1  A2) /< )  *  C  H  MM  SC  A)  -L2>  +  <  AP*X/  b  1 ) 

RETURN 
n  4=  FN  F  (  R 1 ) 

L"T  n3=-D4/C 

M  l  =  -(  Pit  AR)  *22 -GP/C  F?r  <  1+6?) ) 

■v.2=  (  03/  P)  +  (  (  0 lx  D2>  /  P)  +  A2-<  D4 

M3  =(  R?*Dl/p}  + A2S--D2 
M4=C(R2«=D1/R)  +  C  (  A2+  DtD?)  )/R? 

N5=  42  *.  p  1 1  {  A 2-  1 )  K  H  +■  C'l  S\  A)  -  (  0  3/  R  2  >  3) 

ax 
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APPEHDIX  B  (continued) 


0 1 03 0  LFT  M6=<  R1*X  l-X)/(  Pit  2) 

01090  =v«6*A2 

01 100  LFT  K 2—  vi  1  --e €  v)2-(M3-!c'<i4)  J+'lS+'lf 
OHIO  RFTINN) 

01-120  OFF  frgACRl)=S0RCFlt?+4*Dt?-^!'Diel::l*C'?S«^>> 
01130  OFF  FV  BC  R>  =  A  3  1 F  -  3*  «2*  <  '•'  ?/  rO  ♦  A  4 
01  l/iO  OFF  FVCC  P 1)  s  C  1"  CC  FI- 2*D*C3  -SC  A)  )  /F?>  >  /C 
01150  D?57  F\!  DC  ft 1  >  =  <  R 1  -  2  *  T'*  C’5*  S\  A)  >  /  ft? 

OllfiO  OFF  F'lc<  F  O  =  C  P2-  C  ft  1- 2*  P*C3  Sv  A)  )  *  D2>  /  R2t  2 
01170  wrv|  RFAP  O’JMFTAMTFt  3>A2»A3*A4 
0113 0  r^-T*  2005  0*  1  •  09*  •  07*  - •  14 
01190  RFM  t<~ fcp  CHARPF  FFIOHT*  CHAftOF  PFPT1-). 
01200  DATA  40000*  200 
01210  RFM  ftFAR  PAC-F 
01220  DATA  407*220 
09  000  FS1  n 


EXAMPLE  OP  OUTPUT: 


5TV232  14.  54.  19.  04/20/7? 


•/=  40000  DFtJTH=  20  0 

I'/|AP6'  OH  A  |-P“  V  "!  0-3T  11.7255 

P  f:  L  pr  -  ppo  3  pAf7*  407  ft  M~ P=- 5 4, 99 32  ft'F  P  MlCAUs- l 3A3» 31 
FCT)=  A. 283 13  T/ 1H" TA=- 5 . 77373 

r\j  p  1 ;r  riAJA  AT  270 
14.5  4.23.  PASIC  FX=,0"TnM  "PI  Tft 
CH'  A-^IpT 
♦  ♦PB-ADY. 


#3 
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